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Abstract. We construct 9-parameter and 13-parameter dynamical systems 
of the plane which map bi-quadratic curves to other bi-quadratic curves and 
return to the original curve after two iterations. These generalize the QRT 
maps which map each such curve to itself. The new families of maps include 
those that were found as reductions of integrable lattices by 

1. Introduction 

The QRT maps were proposed by Quispel, Roberts and Thompson [H] about 
twenty years ago and provided a fruitful starting point for the held of discrete 
integrable systems. In particular, the de-autonomization of these maps through 
the singularity confinement method led to the discrete Painleve equations, which 
turned out to have very rich connections with other fields, including random matrix 
theory. 

However, although the QRT mappings are considered as the most general family 
of Liouville integrable bi-rational maps of the plane, there exist examples of inte- 
grable maps that do not fall in this class. We consider a generalization of QRT maps, 
which have been introduced under the name of HKY mappings. Their non-QRT 
nature arises from the fact that they preserve a bi-quartic rational integral instead 
of a bi-quadratic one. Actually, the iterates of these systems alternate between 2 bi- 
rationally equivalent elliptic curves, hence their non-QRT nature. Several examples 
of such systems were exhibited in [7] by Hirota, Kimura and Yahagi, where they 
presented third-order mappings, which can be integrated to second-order with bi- 
quartic invariants, hence the initials HKY under which these examples are known. 
Nevertheless, the first researchers that gave examples of non-QRT type mappings 
were Haggar, Byrnes, Quispel and Capel in [6]. For this reason, in this paper, we 
describe these systems as non-QRT mappings. More examples of non-QRT systems 
were generated by autonomisations of discrete Painleve equations [12] , reductions 
on integrable lattices [IT], using an elliptic function solution [18] . solutions of Adler's 
lattice equation [Jj , or by a direct approach [TBI H] ■ 

In this paper, we extend these known examples by providing multi-parameter 
families of non-QRT mappings. Our main idea is to construct birational mappings <fi 
which transforms an integral I(x, y) of a QRT mapping to an involutive homography 
of itself. The construction then leads to mappings that preserve a bi-quartic integral 
and hence lead to non-QRT mappings. This construction is explained in further 
detail in Section 2. In Sections 3 and 4, multiparameter families of examples of non- 
QRT Liouville integrable mappings of the plane are presented. Finally, in Section 
5, we discuss the integration of these systems, followed by conclusions in Section 6. 
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2. From QRT to non-QRT Mappings 
The QRT mapping is denned through a biquadratic invariant 



I{x,y) 



X T A Y 



where 



X T AiY ' 

where X, Y are vectors X = (x 2 , x, 1) T , Y = (y 2 , y, 1) T and Aq, A\ are two 3x3 
matrices, 

§i <, Ci 

Xi [1% 

If these matrices are symmetric, i.e., Ai = Af , the QRT mapping is called sym- 
metric since its invariant is symmetric under the interchange of x and y. If they 
are antisymmetric, i.e., A; = —Af , the QRT is called antisymmetric, otherwise it 
is called asymmetric. In the antisymmetric and asymmetric cases, the QRT is the 
composition i 2 oi% of the non-commuting involutions i\, i 2l where the latter are de- 
fined by the solution of the equations I(x, y) — I(x, y) — and I(x,y) — I(x, y) = 
respectively. So we have 

_ h(y) ~ h(y)x ( x = x 

k ■ { f2{y)-h{y)x , i 2 : < _^ giQc) -92{x)y , (1) 

y = y { 92(x)-g 3 (x)y 

(/i(y),/ 2 (2/),/ 3 (y)) T -(A Y)x(A 1 Y), 

( 9l (x),g 2 (x),g 3 (x)) T = (AtX.) x (AfX). 

In the symmetric case, the QRT is defined by the composition of i\ or i 2 with 
the involution j : x — y, y = x, that clearly preserves I(x, y) since the latter is 
preserved under the interchange of its arguments. So the symmetric QRT can be 
written as the 3— point map: 

__ fijy) - h{y)v 
V h{y)-h{y)y 

The results of ill! suggest a direct approach to finding non-QRT mappings from 
the QRT integral. Specifically, it was shown that if one can find the appropri- 
ate I(x, y) or equivalently the appropriate Ai such that a birational mapping <fi 
transforms I(x, y) to an involutive homography of itself, namely if 

ai — b . . 

Iocj>=- , 2 

cl — a 

then <f> will preserve the bi-quartic / H — and will be an non-QRT map. Note 

cl — a 

that in the case when <f> is an involution, or more generally a periodic map, one can 
get nontrivial non-QRT systems by composing <j) with the QRT involution i\ or i 2 . 
Note that in this case, the non-commutativity of i\ or i 2 with cf> should be checked. 

The authors of [11] identified two cases, which we call Type I and Type II non- 
QRT maps, where 

Type I : I o <fi = —I 
Type II : Io = l/I 
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Equation ((2]) can be mapped to each of these two types. If c = 0, type I can be 
found after substituting I I — b/a. If c ^ 0, we get type II after substituting 
I — » lyj a 2 /c 2 + b/c + a/c. Therefore, the conserved quantities of <f> are I 2 and 
I + 1/ 1, respectively. These two families of mappings at a first glance appear to 
be different, but are related [2J, since if <f> gives a Type / invariant, after defining 

K = ^ — j , we get K o <j) = 1/K. So, from now on, we deal only with Type I 

non-QRT mappings. 

Note also that for both types of non-QRT maps, their even (odd) iteration is a 
QRT mapping, since after two consecutive actions of the map, the bi-quadratic is 
mapped to itself, / — > —I — > I or / — > l/I —> I. In this sense non-QRT maps can 
be described as a "square root" of the associated QRT map. 



3. Two Families of Antisymmetric non-QRT Maps 

Many examples of 3-point non-QRT maps (<pi) that appear in the literature are 
defined for symmetric I(x,y). That is, 4>i = j ° hi, where 

v ,yj , ] : I _ y , and hf = f = id 

y=y \ v=x 

with j of QRT type and hi of non-QRT type (without loss of generality of type I). 

A natural question then arises: for which bi-quadratic I(x, y) is the involution 
j of non-QRT type? If we can find such an I then by composing j with a QRT 
involution i\ or ii of section (2), then we will clearly have an non-QRT mapping. 
The question above is answered in Proposition 13.11 followed by Corollarv l3.H 

In Corollary 13 . 21 the same question is answered for the period 4 map j: x — 
—y,y — x, where again by composing j with the QRT involution i\, we have 
another family of non-QRT maps. 



Proposition 3.1. For 



N(x,y) X T A Y 



when A\ = Aj and Aq — — eAQ , there is I(y, x) + el(x, y) — 0, where e = ±1. 
Proof. When A\ = A\ there is: 

D(y, x)-D(x, y) = Y^iX-X^Y = Y T A!X-Y T ^X = Y T (A^ -Ai)X = 0. 
Also for A = -cAq, 

N(y, x)+eN(x, y) = Y T A Q X+(3, T A a Y = Y T A X+eY T ^X = Y T (^+e^ )X = 0. 
So 

N(y,x) N{x,y) N(y, x)D{x, y) + eN(x, y)D(y, x) 

I{y,x)+el{x,y) = — r + r = ^7 77^7 ; = 

D{y,x) D(x,y) D{x,y)D(y,x) 

D{x,y){N(y,x) + eN(x iy )) 
D(x,y) 2 

□ 
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Corollary 3.1. For A\ — Af and Aq — — eAg, the map 
x — y, 

_ h(y)-f2(y)y , A . 

MV)-MV)* ° r y=My)-f 3 (y)y> (4) 

h{y) - h{y)x 

where 

(hJiJsf = (A a Y)x(A 1 Y), 
for e = —1 preserves I{x,y) |3J) so is a symmetric QRT and for e = 1 is measure- 
preserving with density m(x,y) — (X r ^4oY) and preserves I(x,y) 2 so is a non- 
QRT Liouville integrable map [17j . 

Remark 3.1. Non-QRT mappings presented in subsection (4.2) of [5] and in sec- 
tion (2) of [5] are subcases of the mapping of the corollary (3.1]) . 

Remark 3.2. QRT mappings with A\ = Aj and Ao = — Aq , together with their 
singularity patterns and their connection to discrete Painleve, were studied in [19] 
under the name of antisymmetric QRTs. From the corollary above, it is clear that 
the even(odd) iterations of the mappings of the corollary (3.1\) . with e = 1, are the 
antisymmetric QRTs. In other words, we have shown that the "square root" of the 
antisymmetric QRTs are non-QRT mappings. 

Example 3.1. When: 



/(e-l)« 1 
A =\ -e (e-l)A 

V " £ 7 (e 



















I 





• 


- 1)0 / 












according to corollary (|3.ip there is 

I(x,y) = x 2 y + 7X - e(xy 2 + jy) + (e - l)(KX 2 y 2 + Xxy + fj,), 



, f^n + (e - l)Ay n - 7 



or 



_ «&n + ( e - l) Aa; n ~ 7 
T" X n — i - 



(e - l)«a;2 + 

When e = — 1 mapping preserves the symmetric integral I{x, y) and is a member 
of the symmetric QRT. 

When e = 1, I(x,y) is antisymmetric under the interchange of its arguments 
and the mapping preserves I(x,y) 2 . Hence it is a non-QRT mapping. Moreover 
the even (odd) iterations of <f> preserves the antisymmetric I(x,y). Hence it is a 
member of the antisymmetric QRT 19 . To see the latter, when n is even, note 
that <j> reads: 

X2m+1 = V2m, 2/2m+l = ~X2m + 2/2m - j/VZm, ( 5 ) 



1 According to |15| . a mapping of the plane <f> : (x, y) — > (x, y) is (anti) measure-preserving with 
density m(x, y), if the Jacobian determinant (J) of <f> can be written as J = ^j|^=y. Anti measure- 
preservation corresponds to a measure-preserving and orientation-reversing mapping <j>. For all 
mappings presented in this work, (anti) measure-preservation is proven by direct computation. 
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When n is odd, we have: 

X2m = 2/2m-l, Dim = -»2m-l + Vim-X ~ l/Uim-l, (6) 

Substituting ([6]) into ([5]) we find that the odd x and y satisfy the antisymmetric 
QRT map: 

X2rn+X = ~ %2m-X + Uim-X ~ 1 /V2m-X, 2/2m+l = ~ J/2m-l + ^2m+l - 1 / X2m+X 

Similarly, the even x and y also satisfy the antisymmetric QRT: 

X2m+2 = -X2m + Vim ~ l/Vim, U2m+2 = -J/2m + %2m+2 - l/x2m+2- 

Asking the same questions as above for the period-4 mapping j: x — —y,y = x, 
we find the following corollarjQ 

Corollary 3.2. The mapping 

x = -y, 

- = fi{y) - h{v)x = (iy{5y 2 + C) + ag(2 Sy 2 + 1V A + Q ( 7 ) 

V f2(y)-Mv)x pxy{8 + 1 y 2 )-a{25y 2 + 1 y i + Q 

or 

_ = Py(Sy 2 + 0-ay(2Sy 2 + 7 y 4 + Q 
V -/%y(<5 + 7 t/ 2 )-a(2<V + 7?/ 4 + C)' 
preserves I 2 where I is given by 

X T A Y ( a \ 

i( x >v) = x_ t a y ' with ^° = ^ and 

1 y —a I 

so is a non-QRT mapping. Mapping is anti-measure preserving with density 
m(x,y) — x' j ' A Y ' h ence a Liouville integrable map. 

4. Asymmetric and Symmetric non-QRT Maps 

Here we present a family of non-QRT mappings that corresponds to asymmetric 
QRT case, and then we specialise this result to the ones that correspond to the 
symmetric QRT case. 

Proposition 4.1. For 

I (x,y) = , where A = I 6 e C , M = 

y k a fi ) \ " i . m / m / 

(8) 
with 

a.\ — 2a 2 f — ca5 + (35b — bae /3\ = 2 fa/3 + 67$ — cat — baC, 71 = 2faj — caC, 

Si = 2ma 2 — 2can + 2b/3n — 2ba\ ei = 2ma[3 + 2&7K — 2caX — 2ba[i Ci = 2maj — 2ca\i 
K\ = ma5 — 2 fan + ben — b5X Xi = mae + 6£k — 2faX — bS/i [i\ = ma( — 2/a/x 

there is I o i 2 — I , and I o hi = —I, i = 1,2, where i 2 the QRT involution and 
( __ (<fy 2 + + S)x + 2(k?/ 2 + Xy + /x) 
ftl : I X " 2{ay 2 + j3y + j)x + Sy 2 + ey + 5 ' 




Q 1 


ft 


7i 




ei 


Ci 


«1 


Ai 


Mi 



2 The proof is exactly analogous to the proof of corollary 1 13.11 ) 
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(Sxy 2 + ay + Si)x + 2(n 1 y 2 + X x y + m) 

X = 

2(aiy 2 + f3iy + ji)x + 5iy 2 + eiy + 61 

y = y, 

the non-QRT involutions. 

Proof. Lets consider the function: 

j., _ xx(ay 2 + fly + 7) + x(py 2 + ay + r) + x(<5y 2 + ey + £) + ny 2 + Xy + fx 

xx(ay 2 + by + c) + (x + x)(dy 2 + ey + /) + fcy 2 + ly + m 

(9) 

where a, /?,..., Z, m are parameters and a; = f(x,y) is a function of to be 
determined to be an involution, i.e., 2; = x. Now the solution of the equation 
Iix, y, x) + I(x, y,x) — 0, for x will determine the form of the latter. Explicitly x 
reads: 

_ _ x (y 2 {5 + p) + y(e + a) + C + r) + 2( K y 2 + Ay + fi) 
X ~ 2x(ay 2 + 0y + -f) + y 2 (S + p) +y(e + a) + ( + t 
Substituting this solution back to @, we get / = I(x,y) which is quadratic in x 
and of higher degree than quadratic in y. Clearly, the equation I(x, y) + I(x, y) — 
by construction has as a solution the previously defined (x). A choice of the 
parameters of ([9]), that will lead to biquadratic I(x,y) is p = 6, er = e, r = 0, a — 
d = k = 0, e = M 5 1 = Then and after the scaling (S, e) (5/2, e/2), the integral 
X. T A Y 

I(x, y) = T gets the form above and the solution of I(x, y) + I(x, y) = is 
exactly h\,h2- 

□ 

Corollary 4.1. TTie mapping 

(Sy 2 + ey + S)x + 2(ny 2 + Xy + u) 

x == 

2(ay 2 + fly + j)x + dy 2 + ey + S 

- = 9i{x) - yg2(x) 

92{x)-yg 3 (x) 

and the dual mapping 

(<W + eiy + Si)x + 2(Kiy 2 + Xxy + m) 

x == - 

2(aiy 2 + fliy + 71 )x + c5iy 2 + e x y + §x 

_ = 9i{x) - yg2{x) 

92{x) - yg 3 (x) 

preserve I 2 where I is given by {![), so they are non-QRT maps. Mapping U0\) 
is measure-preserving, whereas mapping is anti measure-preserving, both with 
density m(x,y) — xr ^ Y , hence Liouville integrable maps. 

If we now demand the matrices of ([5]) to be symmetric, we have the following 
corollary: 

Corollary 4.2. The mapping 

x = y, 

2~fy 2 + 25y + 2p + x(fly 2 + ey + S) ( 12 ) 



(10) 



(11) 



y = 



fly 2 + ey + S + 2x(ay 2 + fly + 7) 
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or 

2 ly 2 + 25y + 2^ + y{fly 2 +ey + 5) 

y = = ^ 

fly 2 + ey + S + 2y(ay 2 + fly + 7) ' 

and i/ie duaZ mapping 
x = y, 

2 liy 2 + 2S lV + 2m + x{fl iy 2 + e x y + St) 

y = 

fliy 2 + e iy + <5i + 2x(a 1 y 2 + fliy + 71) 

or 

2 lx y 2 + 25iy + 2m + y{fliy 2 + ay + <5i) 
fty 2 + eiy + <5i + 2y(aiy 2 + fliy + 71) 
preserve I 2 where I is given by 

X T A Y 

I i x > y) = xTA Y i w/iere A 



(13) 











• 


e 








5 


M / 







ft 


7i 


ft 


ei 




71 




Ml 



(14) 



ai = ae - /3 2 ft = 2a<5 - 2/3 7 7l = /3£ - 7 e 

a = ietfj, - 47 2 <5i = 2ftu - 2^/8 m = e M - $ 2 
so they are non-QRT maps. Mapping ilty) is anti measure-preserving, whereas 
mapping U3]) is measure-preserving, both with density m(x,y) — x t^ qY ' hence 
Liouville integrable maps. 

Remark 4.1. For a = 0, fl — 0, 7 = 1/2, S = 0, e = a, yu = 6/2, i/ie integral ( fLfl ) 
iafces i/ie form: 

_ 2ayy + y 2 + y 2 + b 

[V -' V) - 2yy + a(y 2 + y 2 )-ab 

that the associate mapping 112^1 is the non-QRT mapping associated with Q3 found 
mHI]. 

Remark 4.2. For a = k 2 sn(a)sn(b)sn(a — 6), fl = 0, 7 = —sn(a — b), 6 = 0, 
e = 2(sn(a) — sn(b)), // = sn(a)sn(b)sn(a — b), the integral |i^| ) tafces i/ie form: 

1% v) = 

((1 + k 2 y 2 y 2 )sn(a)sn(b) - y 2 - y 2 )sn(a - b) + 2yy(sn(a) - sn(b)) 
((1 + k 2 y 2 y 2 )sn(a)sn(b) + y 2 + y 2 )(sn(a) — sn(b)) + 2yysn(a — b)(k 2 sn(a) 2 sn(b) 2 — 1) 

i/iai the associate mapping ilty) is the non-QRT mapping associated with QA found 
in \U\. 

Example 4.1. For a = 1, fl = S = 0, I given by (jT4j) takes the form: 



, , _ x 2 y 2 + 1 {x 2 + y 2 ) + /i _ j/V+7(j/ 2 + y 2 ) + ^ 

i( x iy) — i or i\y>y) — > 

- yy 

and from (1121) there is: 



?yy = — 5— — (15) 
y 2 + 7 

that preserves I(y,y) 2 . Inspired by [8], a discrete Lax pair for (fTSf is given by 

L„*„ = A*„, *„+i = M„*„, (16) 
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where 




/ V / a I «■ — ± if t * I 

L n (A) = [ ^ n _ 2 n _ 2 r_l\n V+liVn-i+vl.) 



Vn-lVn 



M„(A)=I \ a „2 / 7 n n_. M+7^+^-i(7+l/S) )■ 

The equation ([15]) arises as the compatibility condition L„ +1 M„ — M„L„ = for 
the linear system (Q2 



4.1. Integration of the non-QRT mappings. In exactly the same manner as 
in [13] . and [9], by integration of a non-QRT mapping we will mean the parametri- 
sation of the corresponding invariant curves. 

In this setting, the integration of the symmetric QRT, according to Veselov 
[17] was first considered by Euler, but of course not inside the modern context of 
integrability. A more transparent and explicit method was presented by Baxter [3] . 
Recently, the authors of [T3], [5], independently presented the integration of the 
asymetric QRT. 

To proceed with the integration of the non-QRT mappings Q and (J7J), we should 
first make the following observation. Let us consider 

77 ■ 1 — ( -nn ^-"-!^ ^" 

■l\Xn—l,X n ,rl) — { I) T , 

then 

I{x n , x n+ i;n + 1) - HXn-i^Xn] n) = (-1) + h ^ — 

\A n AiA n+ i A^^iAn 

that is satisfied by mappings and {?]). Then the invariant curves that the HKY 
under consideration lies is: 

XLi^oX,, - (-lJ-A-^xAiXn = 0, (17) 

where K is the integration constant associated with the initial conditions (x—i, xq), 
K = I(x-i,xq). Then, after defining M = Aq — (— l) n KAi, the way to parametrise 
X^ , _ 1 MX TS = 0, is exactly the one presented in [13] for the asymmetric QRT. 

5. Conclusions 

In the previous sections, multiparameter mappings of non-QRT type were pre- 
sented. By construction they preserve a bi-quartic invariant that can be consid- 
ered as the square of a QRT bi-quadratic one. Apart from the preservation of a 
bi-quartic integral, these systems turn to be (anti) measure-preserving and hence 
Liouville integrable. 

Commenting on the search for general solution of such integrable systems, it 
should be noted that by fixing an initial condition say (x_i,Xo), we have K = 
I(x-x,Xq) and equation (|17p defines, in general, two elliptic curves: one for n odd 
and the other for n even. The non-QRT system at each step maps one curve to 
the other and vice versa and since the latter is bi-rational, these two curves are 
birationally equivalent and hence they should have the same j-invariant. In this 
sense we should not expect to find the general solution of a non-QRT map as an 
addition formula on one elliptic curve since the latter simply alternates between 
two birationally equivalent elliptic curves. Also in this setting by construction, by 
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finding a non-QRT map we have two birationally equivalent curves and the non- 
QRT map is the birational relation between these two curves. Note also that these 
are special classes of birationally equivalent curves since the birational relation that 
maps one to the other and the inverse of it are the same birational map and the 
one we call the non-QRT map. 

It should be pointed out, that although general multiparameter families of these 
systems are presented, we have by no means exhausted all possble non-QRT map- 
pings of the plane. A natural extension of this work is to consider systems whose 
higher iterates are QRT maps. If the third iteration is a QRT mapping, then they 
preserve I 3 , with I associated to the QRT bi-quadratic invariant. This will be 
considered elsewhere. 
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